THE PARABOLIC TWO-PHASE MEMBRANE PROBLEM: 
REGULARITY IN HIGHER DIMENSIONS 
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Abstract. For the parabolic obstaclc-problcm-likc equation 

Au - dtu = X+X{u>o} - ^-X{u<o} , 
where A+ and A_ are positive Lipschitz functions, we prove in arbitrary finite 
dimension that the free boundary d{u > 0} U d{u < 0} is in a neighborhood of 
each "branch point" the union of two Lipschitz graphs that are continuously 
differcntiable with respect to the space variables. The result extends the elliptic 
paper 1101 to the parabolic case. The result is optimal in the sense that the 
graphs are in general not better than Lipschitz, as shown by a counter-example. 
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Figure 1 . Example of a Stationary Branch Point 



1. Introduction 

1.1. Background and main result. In this paper we study the regularity of the 
parabolic obstacle-problcm-like equation 



where T < +oo, A_(- > 0, A_ > are Lipschitz functions and C R" is a given 
domain. The problem arises as limiting case in the model of temperature control 
through the interior described in [U 2.3.2] as hi, h2 0. 

We are interested in the regularity of the free boundary d{u > 0} U d{u < 0}. 
As the one-phase case (i.e. the case of a non-negative or non-positive solution) is 
covered by classical results, and regularity of the set {u = 0} n {Vu ^ 0} can be 
obtained via the implicit function theorem (see Section ?? for higher regularity), 
the research focusses on the study of d{u > 0} fl d{u < 0} n {Vu = 0}. 
In the stationary case — the two-phase membrane problem — the authors proved 
f [TT] and [TU]) that the free boundary d{u > 0} U d{u < 0} is in a neighborhood of 
each branch point, i.e. a point in the set H, n d{u > 0} n d{u < 0} n {Vu = 0}, the 
union of (at most) two C^-graphs. Note that the definition of "branch point" does 
not necessarily imply a bifurcation as that in Figure [TJ 
We formulate the main result in this paper. 

Theorem 1.1. Suppose that 



(1.1) 



Au - dtu = X+X{u>o} - ^-X{«<o} 



in (0,r) X n, 



< Ainin < inf min(A+, A_), sup max(|VA+|, |VA_|, l^tA+l, |atA_|) < +oo 

QUO) Qi(0) 
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and that u is a weak solution of 

Au ~ dtu A+X{«>o} - A_X{„<o} in (3i(0) ; 

here Qi(0) is the parabolic cylinder (—1, 1) x i?i(0). 
Then there are constants cr > and rg > such that 

(1.2) w(0) = 0, |Vw(0)| < (T, pardist(0,{w > 0}) < cr and pardist(0, {u < 0}) < cr 

imply d{u > 0} n Qro(O) and d{u < 0} n (5ro(0) being graphs of Lipschitz func- 
tions (in some space direction) that are continuously differentiate with respect 
to the space variables. The constants cr, rp , the Lipschitz norms and the modu- 
lus of continuity of the spatial normal vectors to these surfaces depend only on 
infg^jo) niiii(A+, A_), the Lipschitz norms of X±, the supremum norm of u and the 
space dimension n. 

Moreover the regularity above is optimal in the sense that the graphs are in general 
not better than Lipschitz. 

Corollary 1.2. Suppose that 

< A,„in < inf min(A+,A_), sup max(|VA+|, |VA_|, |9tA+|, |atA_|) < +oo 

Qi(o) 

and that u is a weak solution of 

Au - dtU ^ \+X{u>o} - A_X{m<o} inQi{0). 

Then there is a constant ro > such that if the origin is a branch point, then 
d{u > 0} n Qra{0) and d{u < 0} H (5,o(0) f^^G graphs of Lipschitz functions (in 
some space direction) that are continuously differentiable with respect to the space 
variables. The constant ro, the Lipschitz norms and the modulus of continuity of 
the spatial normal vectors to these surfaces depend only on infQj(Q) inin(A-(-, A_), 
the Lipschitz norms of X±, the supremum norm of u and the space dimension n. 

As to the proof we extend the method of [10] to the parabolic case. There is 
however a difficulty as the time derivative dtU is in general not continuous, so that 
it is not possible to apply directly the comparison principle. We deal with that 
problem by a two-stage proof of directional monotonicity. 

2. NOTATION 

Throughout this article R" will be equipped with the Euclidean inner product 
X ■ y and the induced norm |a;| , Br{x^) will denote the open n-dimensional ball of 
center , radius r and volume r" w„ , B[, (0) the open n — 1-dimcnsional ball of 
center and radius r , and the i-th unit vector in R" . We define Qr{t^ :X^) := 
(t° — , <° + ) X Br (x" ) to be the cylinder of radius r and height 2r^ , Q~{t^ ,x^) : = 
(t° - r^, t°) X Br[x°) its "negative part" and x^) := i° + r^) x Br{x°) its 

"positive part". When omitted, a;° (or (t°,a;°), respectively) is assumed to be the 
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origin. Moreover let 9parQr(<°, a;°) := -r'^ ,t° + r'^)xdBr{x°)\J{f -r"^} y. Br{x°) 
denote the parabolic boundary of Qr{t'^,x'^). Let us also introduce the parabolic 
distance pardist((i, x), A) := inf(,,j,)g^ y/\x — yp + |i — s| . Given a set A C R"'*""'" , 
we denote its interior by A° and its characteristic function by XA ■ By Vu we 
mean the gradient with respect to the space variables. In the text wc use the n- 
dimensional Lebesgue- measure and the m-dimensional Hausdorff measure 7i™. 
Finally, C'^'^ :— H^'^ denotes the parabolic Holder-space as defined in [7]. 

3. A SUPREMUM-MEAN-VALUE ESTIMATE 

In this section we show that at branch points the time derivative dtu, in general 
a discontinuous function, satisfies a sup-mean- value estimate. 

Lemma 3.1. Let Q~^(t°,x°) C (0,r) x n and let A+,A_ be non-negative and 
Lipschitz continuous with respect to the time variable. Then each solution u of 
hl.l]] satisfies 




sup \dtur,\ = sup \dtu\ <C\r'+\ r-"-^ / \dtu\ 

Qi Q^(t«.x°) y V JQ-^{t°.xO) 

Proof. Using the scaling invariance of the equation with respect to the scaling 

Ur{t, x) = r~'^u{t° + r'^t, x° 4- rx) 

we may assume that r = 1/2, t° = and = 0. 
Let H{t,x,z) = X+{t,x)x{z>o} - ^-{t,x)x{z<o}- For 

u \ P,r u \ u{t + T,x) -u{t) 
v(t, X) := u[t, x) :— 

T 

and ?7 G ^^((-1, 1); W^''^{Bi)) such that r/ = on (-1,0) x dBi, we calculate 
= -f_J^^i^dlH{t,xMt,x)), .se(-l,0). 

Here 

dlH{t, X, u{t, x)) = A+(t" + rh, x° + rx)dlx{u>o} - A_(t" + r\ x° + rx)dlx{u<o} 

+X{u{t°+r^{t+T),x°+rx)>Q}dt\+ — X{u{t° {t+r) ,xO +rx)<a}dl \- ■ 

Testing with r){t,x) := C{x)(l)^{t)xasx{v{t,x)-k,Q) where fc > 0, C e Cl'^{Bi) and 
e C"'H-1, 1) such that (j){t) e [0, 1] and 

m ■■= 

and observing that 

max{v{t, x) — k, Q)dlH(t, x, u{t, x)) > —Cir^ max(z;(i, x) — k, 0) 
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we obtain 

(3.2) sup / 0^(.s)C^max(w(s,-)-fc,O)^+ / / (l>^(:^\\/ max{v ~ k,0)f 



<C2 / [max(w-fc,O)^(r|VCr| + 0|at0|C)+HrC max(i;-fc,0)] 

J-l JBi 

From the proof of [3 Theorem 4.7] we infer that 



(3.3) supv <C3 + { I 

Ql/2 

Testing with ri{t, x) := C,"^ {x)4)'^ [t) max(— i'(t, x) — fc, 0) where > 0, we obtain in a 
similar way that 

(3.4) sup(-w) <C3lr^+ [ vA . 

Letting t ^ and scahng back wc obtain the statement. □ 

4. NON-DEGENERACY AND REGULARITY OF THE SOLUTION 

Lemma 4.1 (Non-Degeneracy). For every Q2r{t''\x'-*) C (0,r) x il the following 
holds: 

l)lf(t°,x°)£d{u>0}, then sup u> — inf A+ . 

2) If {t° , x°) e d{u < 0} , then inf u< — - inf A_ . 

Proof. We choose a sequence {u > 0} 3 (t™, x™) — > {t°, x°) as m — > oo . Supposing 
that supg-^jm ^m-j u < ^ iuf (jo ^0 ) A^. , the comparison principle yields that 
u{t,x) < v{t,x) := {^-f^ + -^\x - Xj-n\^) infQ^(to_^o) A+ in (37(t™,a;™) , a contra- 
diction to the fact that w(i™,x™) > . 

The estimate for infg-j-jo ^.o) u is obtained the same way, replacing u by —u and \+ 
by A- . ' □ 

Lemma 4.2. Let A+, A_ £ Cio'c ((0' ^) ^ Then each solution u of (Op satisfies 
the following: 

2)dtVuGLl^{iO,T)xn). 
Proof. 1) follows from Lemma l3. II 

2) follows from ()3.2p with = and from the analogous estimate for max(— w,0). 

□ 
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Corollary 4.3. For every Q2r{t°,x^) C (0,T) x O, there exists a constant co > 
depending only on n and ||9t?i||^oo(Q^(-to^o-|) such that 

u > in {tP ,x^) implies u > in Qcar{t'^ ,x'^) i O'lT-d 

u < in Q~ {t'^ ,x'^) implies u < in Qcori't'^ : x'^) . 

Proof. Suppose towards a contradiction that u{t^,x^) < for some {t^,x^) G 
Then there is a point (^^a;2) € d{u < 0} n Qtor{t°,x^). Apply- 



ing Lemma |4?T] at (t , a; ) with respect to the cyHnder Q(i-co)r{t ,x ) yields a 
contradiction to Lemma 14.21 1) provided that cq has been chosen small enough. 
The second estimate is proved in the same fashion. □ 

Proposition 4.4. Let A+,A_ € (7[°'^((0,T) x Q). Then each solution u of il.l]) 
satisfies Vu € C^^^'^ {{0,T) x fl), that is, the gradient is Lipschitz continuous with 
respect to the space variables and Holder continuous with exponent 1/2 with respect 
to the time variable. 

Proof. Let us first show that for any e G dBi, (A — 9t)(max(9eU, 0)) > — C and 
(A — 9t)(max(— 9eM, 0)) > — C in il. We give a formal proof that can be made 
rigorous translating everything into a weak formulation. In {deU > 0}, 

deU 



-(A+H"-i [({Vu ^ 0} n d{u > 0}) + A_H"-i L({Vu ^ 0} n d{u < o})) 



\\/u\ 

+(?eA+X{u>0} ~ de^-X{u<0} > . 

As deU is continuous, wc obtain (A — 9t)(max(9eU, 0)) > — C. 
Considering — e instead of e wc obtain also (A — 9t)(max(— 9,.^, 0)) > — C. But 
then the "almost monotonicity formula" Theorem I of [5] applies and we proceed 
as follows: at each point x") € {u ^ 0} fl {Vm = 0}, we obtain from the almost 
monotonicity formula that "S/deU is bounded at {t'^,x^) by a locally uniform con- 
stant. 

At each point {t'^,x'^) E {u ^ 0} D {Vu ^ 0}, we obtain in a similar way that 
for every e_LVu(f'^, x*^), \\'deu{t'^ , x'^)\ is bounded by a locally imiform constant. 

^1 = |vu(fV»)| ■ "^^^^ -diiu{t°,x°) = -A+X{„(to,2;0)>o} + \-X{u{to.xO)<o} - 
dtu{t^, a;") + X]j=2 djju{t^, x^) is by Lemma |4?2] bounded by a locally uniform con- 
stant. □ 

Corollary 4.5. C'+^{d{u > 0} U d{u < 0}) = 

Proof. First, we obtain from Lemma H?T1 Lemma and Proposition l4.4l that there 
exists a locally uniform constant c > such that for Q2r{s, y) C (0, T) x fi, 

^ > c> If (s, y) e a{« > 0} 
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and Cn+i{Q^) > c > if (s, y) G d{u < 0} . 

Since X{u>o} * XQ./'C"+i(Q^) ^ X{„>o} in Lii„^((0,T) x 1]) as r ^ and the 
analogous fact holds for X{«<o}j we obtain that X{„>o} > c > £"+^-a.e. on 
9{u > 0} and X{u<o} > c > £"+i-a.e. on ^{m < 0}. Thus > 
0} U d{u < 0}) =0. □ 

5. Vanishing time derivative 

As a corollary of Lemma |3. II we obtain now that at points at which the blow-up 
limit depends only on the space variables, the time derivative dfU - in general a 
discontinuous function - attains the limit 0. 

Corollary 5.1. Let Q2r{t^,x^) C (0,r) x Q, and suppose that for a sequence of 
solutions Uk in (0, T) x f2 

Ur^.x) = i-k^'^Ukit^ + rlt,x^ +rkx) u^{x) in L\^^{W+^) as ^ . 



Then 



sup |9tUfc|->0 



as rfc ^ 0. 

Proof. The statement follows from Lemma lOI and the fact that dtUr^. converges to 
in L[^Qp(R"+^) as ru 0. The L^-convergence in turn may be shown as follows: 
as dfUk is by Lemma [3.11 bounded in L°°{Qr{t'^ ,x'^)), it is sufficient to prove a.e. 
convergence. For {s,y) G {mq = 0}*^ we obtain from Lemma 14. II that Uj.^ = in 
Qs{s, y) for some (5 > and large k. For (s, y) e {uq > 0} U {uq < 0}, converges 
in C^{Qs{s,y)) for some 5 > as fc ^ oo. Moreover we know from Corollary 14.51 
that C'+^{d{uo > 0} U d{uo < 0}) = 0. It follows that dtUr^ converges r'+^-a.e. 
to dtUQ. □ 

6. Directional monotonicity 

In a first stage, we show that if the solution is close to the one-dimensional 
solution 

(6.1) h{x) := ^^max(xi,0)' - min(2;i, 0)^ . 

then it is increasing in a cone of spatial directions. Later on we will extend the 
result to a cone of tempo-spatial directions. 

Proposition 6.1. Let < Xmin < infQi(o) min(A_|-, A_), h as in i6.1]) . and let 

e G (0, 1). Then each solution u of in Qi{Q) such that 



distio=((_i^i).vyi,oo(B^))(u, h) <6 : 



48n 
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and 

sup max(|VA+|, |VA_|) < S 
Qi(o) 

satisfies e~^deU — |u| > in Qi/2(0) for every e G dBi{Q) such that ei > e; here 
ei denotes the first component of the vector e. 

Proof. First note that e^^deh - \h\ > in Q2(0). It follows that 

(6.2) e'^deu - \u\ > -36e-^ in Qi(0) 

provided that dist^oo^.i i-j-vi/i oof^^)) (it, h) < 5. Suppose now towards a contradic- 
tion that the statement is not true. Then there exist A+, A_ € (Amin, +00), (t* , x*) £ 
01/2(0), e*, and a solution u of p.ip in Qi(0) such that dist£,=c((_i i^.^yi, 00 h) < 

sup niax(|VA+|, |VA_|) < S, 

Qi{o) 

e* > e and e~^de'u{t* , x*)~\u{t* , x*)\ < 0. For the positive constant c to be defined 
later the functions v :— e'^de'U— \u\ and w := e~^de*u— \u\ + c\x — x*\'^ — c{t~ t*) 
satisfy then the following: in the set D := Qi(0) n {v < 0} n {t < t*}, 

Aw - dtw <2nc + c- A+X{„>o} - A_X{n<o} 
+e-^{X+ + ■ e*W"-i [{{u = 0} n {Vu ^ 0}) 

+£^^{X{u>o}de*K - X{«<0}9e*A_) 

where Vx = As 

i^x-e* <0 on {u = 0} n {w < 0} = {li = 0} n {e~^de*u < 0} , 

we obtain by the definition of S that w is supercaloric in D provided that c has 
been chosen accordingly, say c A,iii,i/(4n). It follows that the negative infimum 
of w is attained on 

5par-D C (aparOl(O) H {t < t*}) U (Qi(0) H a{« < 0}) . 

Consequently it is attained on {t < t*} O 9parQi(0), say at the point {t,x) G {i < 
t*} n 9par(3i(0). Since pardist((f, S), (i*, a;*)) > 1/2, wc obtain that 

e~^de»u{t, x) — \u{t, x)\ = v{t, x) = wit, x) — c\x* — + cit — t*) 

< -c/4 = -A,„i„/(16n) . 
But this contradicts (|6.2p in view of 6 = . □ 
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7. The set of non- vanishing gradient 

In the sequel we are going to need higher regularity of the level set {m = 0} n 
{Vu 7^ 0}. Higher regularity can be obtained in a standard way using the von 
Mises transform: 

Lemma 7.1. The set {m = 0} n {Vu ^ 0} is locally in (0, T) x fl a -surface and 
dtu is continuous on that surface. 

Proof Let {t^,x°) e {u = 0} n {Vu 7^ 0}. We may assume that Vu{t",x^) = 
diu(t'^,x^) and that in Qs{t'^ , x'^), u is strictly increasing in the xi-direction and 
{u = 0} is the graph of a function, say xi = g{t, x') for (i, x) € Qsit^, x°), where g g 
C"{(t° -S^,t" + d^);C\B'g{x°))). It is sufficient to prove that 5 G (Q5/2 , ) ) ■ 
To do so, we use von Mises variables, i.e. 



y = u{t, xi , x') and xi = v(t, y,x') . 



A calculation assures that 



-l-|V'w|2 \ A'v V'vV'dyV dtv 



{dyVf J dyV [dyVf OyV 

^ { X+{t,v{t,y,x'),x'), y>0 
\ -\^{t,v{t,y,x'),x'), y<0 



Thus 

dtv - aij{Vv)di-jV = f{t,y,x')dyV 



-X+{t,v{t,y,x'),x')dyV, y >0 
X-{t,v{t,y,x'),x')dyV, y <0 . 



Provided that 6 has been chosen small enough, |V'w| < 1/2, < dyV < 1/2 and the 
above equation is uniformly parabolic. Moreover, 

dtdl^v - a,, (Vt;)9y - + h, y, x')dkd^v 

OPk 

- fit, y, x')dyd^v + dyv{t + h, y, x')d!^fit, y, x') 

where Zh = 9{t, y, x')Vv{t + h, y, x') + (1 - 0{t, y, x'))Vv{t, y, x') and 9{t, y, x') G 
[0, 1]. Calculating 

div = l/diu, dfV = —dfu/diu, diV = —diu/diu for 2 < i < n, 

dijV + diivdjU + dijvdiU + divdiju + duvdiudjU — for 2 < i, j < n 
diivdiu + diivdiudiU + divdnu ~ 0, duv = —duu/ (diu)'^ 
shows that all spatial second derivatives of v and dtV are bounded. Thus f{t, y, x'), f{t, y, 



and ^2§^d^jv{t + h ,y,x') are bounded uniformly in h, and we obtain from |5] 



that dfV is uniformly Holder continuous with respect to h and that dtU is Holder 
continuous in Q5/2(i°) 2:°). □ 
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8. Global solutions 

In this section we extend our characterization of eUiptic global solutions [21 The- 
orem 4.3] to the parabolic case. We are going to need the following version of the 
Caffarelli-Kenig monotonicity formula of [3]: 



Theorem 8.1. Let 



<i>(r, w) := ^/(r, max(u', 0))/(r, max(— w, 0)) 



where 



I(r,v) / / \Vv\^G{t,x) 

and G is the backwards heat kernel 

G(t,x) = (M-t)r/^exp(^). 

Ifu\.a,yi{w, 0) and max(— ly, 0) are continuous subcaloric functions, then r n- $(r, w) 
is non- decreasing, and $(cr, lu) = <^{p,w) for some < p < a implies that either 

(A) \/ max(w, 0) = m -cr^ < i < or V max(-w, 0) = m -ct^ < < < 0. 
or 

(B) max(w, 0)(9t — A) max(w, 0) = and max(— tu, 0)(9t — A) max(— lu, 0) = m 
— (T^ <t<Oin the sense of measures. 

Proof. For v :~ max(w,0) (or v := max(— ii;,0), respectively) we calculate 



I{r,v)^-]-l I G{t,x){dt-Ay+ f I G{t,x)v{dt-A)v , 

I'{r,v) >2r [ \\/v\^G{~r^,x) . 
Jr" 

In what follows we assume that /(?■, v) ^ 0. It follows that 

I'i^^^) ^^^ /R,.|V^(-r^x)pG(-r^x) 
I{r,v) ~ J^„v'^{—r^,x)G{—r'^,x) 

In the case I'{r,v) ^ the inequality is strict unless J^r2 Jf^n v{dt — A)w 
Consequently <i>(r, w) = 0, or else 



^{r,w) ~ r 



^ ^ ^ J^„\Vma.^iw,0)\^Gi~r^x) 
Jj^„ max(u', 0)^G{—r^,x) 



|Vmax(-u;,0)|2G(-r2,a;) 
max(— w, 0)'^G{—r'^,x) 

where the inequality is strict unless both /°^2 J^in max(ii;, 0){dt — A) max(w, 0) = 
and /°^2 max(— w, 0)(9t — A)max(— w,0) = 0. Moreover, by [31 Corollary 
2.4.6], the right-hand side is non-negative. □ 
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Lemma 8.2. Let be solutions of U.l\) in R"+^ with such that = in 

{t < 0} and v^,v'^ have polynomial growth with respect to the space variables. Then 
v^ = u2 in R"+i. 

Proof. Multiplying the difference of the two equations by {v^—v'^)W where W{t, x) = 
G{t — T, x) and integrating, we obtain for each < T < +oo, < S < T and H 
defined in Lemma 13.11 that = 



S S 

[ W[\\/{v^ ^v^)\^ + {H{v^)-H{v^)){v^ ^v^)] - ^ [ [ {v^-v^)^dtW 

JR" 2 Jq Jji„ 



s 



> 

- 2 



U W{S){v\S)-v'{Syf + U f {v'-vy[-d,W-AW] 

WiS)iv\S)-v'iS)r . 



R" 

□ 



Lemma 8.3. Assume that w is a backward self-similar solution with constant co- 
efficients A+,A_, i.e. 

w{e'^t, Ox) e'^w{t, x) for all e > Q,t < Q and X . 

Then Viy ~ Q on {w ~ 0}. 

Proof. First, the self-similarity implies that 

(8.1) deuiX^t, Xx) = Xdew{t, x) for all e e dBi,\ > 0, t < and .t e R". 

Consequently the function r i— > $(r, dew) of the monotonicity formula Theorem l8.ll 
is constant in (0, +oo), implying by Theorem 18. II that either 

(A) V ma.x{deW, 0) = in {t < 0} or V max(-(9eW, 0) = in < 0}. 
or 

(B) ma.x{deW,0){dt — A) max(9eU', 0) = in {t < 0} and max(— SeW, 0)(9t — 
A) max(— i9eW, 0) = in < 0} in the sense of measures. 

Suppose now towards a contradiction that there is a point (t^ , x^) E {t < 0} {w ~ 
0} n {Vw ^ 0} and denote ^ \lZpP)\ ^nd let Q^{t\x^) such 

that di,ow > in Qi^{t^,x^) and {w = 0} n is a C^-surface. In the case 

i^" • e ^ 0, 

\idt~A)dew\iQ^it\x^))^\\+ + \^\[ [ \e-u\drr-Ut^Q . 

Jt^-K^ J B^{x^)n{w{t)=a} 

Thus (A) holds. From jHH]) we infer that deW > in {t < 0} if e • > and 
deW < in {t < 0} if e • z^" < 0. Hence deW = in {t < 0} for all e±z/°. As in [1 
p. 844] we may write 

w{t,x) = -tf{^^) 
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and calculate the 2-parameter family of solutions of the ODE which /(^) = w{—l, £,) 
satisfies in (0, +00), 

/(e) = A+ + Ci(C2_2) 
+ C2 ^-2Ce«'/4 + ie - 2) e^'/-^ ds^ in {/ > 0} 

and 

fiO = -A- 

+ C4 |^-2ee«'/4 + (e-2)J' e^'l^ dsj in {/ < 0} . 

As w has polynomial growth towards infinity wc conclude that = C2 = C4 and 
that 

/(e) = A++Ci(e'-2) in{/>0} 

and 

/(e) = -A-+C3(e'-2) in{/<0}. 

If /(a) = and /'(a) 7^ for some a £ R then Ci = C3 = -A+/(a2 - 2) 
A_/(a^ — 2), a contradiction. Therefore /(a) = implies /'(a) = and a = 0. It 
follows that Vw = on {w = 0}. □ 

Theorem 8.4. Let w be a global solution with constant coefficients A+,A„ such 
that dtw and D^w are bounded, and suppose that the origin (in time-space) is a 
branch point of w. Then after rotation 

w{t,x) w*{t,x) A+max(.T„,0)V2 - A_max(-a;„, 0) V2 for e R"+i . 
Proof. 

Step 1: Let us first assume that w is a backward self-similar solution. By Lemma 
18.31 Vw = on ui = 0. But then zi := Taax{w,0) and Z2 := max(— WjO) are 
in {t < 0} non- negative backward self-similar solutions. Concerning those, it has 
been shown in Pl Lemma 6.3] and [U Theorem 8.1] that cither Zj is a half-plane 
solution of the form Zj{t,x) = A-|-/2max(2; • e,0)^ for some e € or Zj{t,x) = 

—agt + X]"=i c-i^i with non-negative constants Oi, < i < n. In the latter case the 
symmetry of Zj implies that Zk = in {t < 0} for k ^ j, and by Corollary 14.31 the 
origin cannot be a branch point. 
It follows that after rotation 

w{t,x) = w*{t,x) for i < . 



Step 2: In the case of a general solution w as in the statement of our theorem, 
we consider the blow-up up wq of w at the origin and the blow-down w,^. By 
the non-degeneracy Lemma [4.11 and [T21 Theorem 4.1], both wq and Woo satisfy 
the assumptions of Step 1. Thus both wq and Woc are after rotation of the form 
A4-max(x„, 0)^/2 — A_max(—x„, 0)^/2 for t < 0, and the monotonicity formula [T^] 
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implies that w is backward self-similar. But then it follows from Step 1 that after 
rotation 

x) = w* [t, x) for i < . 
Last, wc apply Lemma [5^ to obtain the same for t > 0. □ 

9. Uniform closeness to h 

We are now ready to prove uniform closeness of the scaled solution to the global 
solution h of (|6.ip , assuming that we are in the setting of Theorem 11.11 

Lemma 9.1. Let u he a solution of hl.l]) in Qi{0). Then, given 5 > Q, there are 
constants rg > 0,(7s > (depending only on inlg^j-Q) min(A+, A_), the Lipschitz 
norms of X± , the supremum norm of u and the space dimension n ) such that the 
following holds: 

If r & (0,r5] , u{s,y) = , |Vw(s,y)| < agr, pardist((s, y), {u > 0}) < asr and 
pardist((s, ?/), {m < 0}) < asr for some {s,y) S Qi/2(0) then in Qr{s,y), the 
solution u{s + + •) is Sr^-close to a rotated version h of the one- dimensional 
solution h defined in i6. 1]) . more precisely 



r ^ sup \u{s + -,y + -) — h\+r ^ sup \\/u{s + ■ ,y + ■) — S/h\+ sup \dtu{s + ■ , y + ■)\ 
QAO) QAo) QAO) 

< S. 

Proof. Suppose towards a contradiction that the statement of the lemma fails. Then 



for some 6 > there exist aj 0,rj 0, {s^,y^) — > (s°,y°) G Q1/2, a sequence 
Uj of solutions such that {s^,y^) G Qi/2(0), Uj{s^,y^) ~ 0, \\7uj{s^ ,y^)\ < ajrj, 
pardist((s^, j/^), {lij > 0}) < crjrj, pardist((s^ , y-'), {uj < 0}) < ajrj and 

rj"^ sup \uj{s^ +r'j-,y^ +rj-)-h{rj-)\ + rj^ sup \Vuj{s^ +ij-,y^ +rj-) -Vh{rj-)\ 
Qi(o) Qi(0) 

+ sup \dtUj{s^ + r'j-,y^ +rj-)\> 5 
Qi(o) 

for all possible rotations h of h. 
We may define 

"-^jy ~r ^- , ' J 



rr , ^ Uj{r]t + s^,rjX + yi) 



and arrive at 

(9.1) \\U,-h\\wi.^iQO>S, 
for all possible rotations h of h. 

Observe that Uj is a solution of (|l.ip in Qi with respect to the scaled coefficients 
X+{r]t + si,r,jx + y^) and A_(r2t + ,rjX + y^). Since C/j(0) = 0, |V[/j(0)| < 
aj, pardist(0, {C/j- > 0}) < CTj, pardist(0, {C/j < 0}) < (jj and the derivatives 
D'^Uj,dtUj are uniformly bounded, we obtain by standard compactness arguments 
a global limit solution Uq of (jl.ip in R" with respect to A+(s", and A_(s°,y°) 
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which satisfies G d{Uo > 0} D d{UQ < 0} n {VC/q = 0}. By Theorem[831 Uo = h 
where h is a rotated version of h. Thus Uj and VUj converge in Qi uniformly to 
h and Wh, respectively, and by Corollarv l5.ll dfUj in L°°{Qi) as j oo. We 
obtain a contradiction to (19.11). □ 



10. Continuity of the time derivative 

Assuming once more that we are in the setting of Theorem II. 1[ we show in the 
present section that the time derivative of the solution is continuous in a suitable 
neighborhood of the origin. 

Proposition 10.1. Let u be a solution of in Qi. Then there are positive 

constants r and a (depending on inf^^ min(A+, A_), the Lipschitz norms of X±, the 
supremum norm of u and the space dimension n ) such that the following holds. If 
u(0) = 0, |Vu(0)| < df, pardist(0, {u > 0}) < af and pardist(0, {u < 0}) < af then 
each blow-up limit at a point (t^,x^) G Qfn{u = 0}n{Vu = 0} is time-independent. 

Proof. Let us consider {t^,x^) E {u = 0} O {Vu = 0}. As the statement of the 
Proposition is by Theorem 18.41 true when {t^,x^) is a branch point, we may as- 
sume that w > in some neighborhood of {t^,x^). From Lemma 19.11 fwith S := 
infQj min(A-)-, A_)/(96n)) and Proposition IG . 1 1 wc know that u is non-decreasing, say 
in the direction e for every e close to x„ in Qf and that \dtu\ < infg- min(A+, A_)/2 
in Qf. 

From 12, , Theorem 4.1] we infer now that each blow-up limit z at {t^ , a;^) is a non- 
negative backward self-similar solution. Concerning those, it has been shown in [51 
Lemma 6.3] and Theorem 8.1] that either z is a half-plane solution of the form 
z{t,x) = A+(i^,a;^)/2max(a;-e,0)^ for some e G dBi, or z{t,x) = —aot + J27=i ^i^i 
with non-negative constants ai,Q < i < n and oq < A+(t^, x^)/2. In the latter case 
the symmetry of z contradicts the fact that z is non-decreasing in every direction 
e as above. Consequently dtz = in {t < 0}, and Lemma [8.21 and CoroUarv 15.11 
imply that dtu{t^,x^) =0. □ 

Corollary 10.2. Let u be a solution of il.l]) in Qi. Then there are positive con- 
stants f and a (depending on inf^^ min(A+, A_), the Lipschitz norms of X±, the 
supremum norm of u and the space dimension n ) such that the following holds. If 
u(0) = , |Vu(0)| < af, pardist(0,{u > 0}) < af and pardist(0,{u < 0}) < af 
then dtu is continuous in Qf. 

Proof. The corollary follows immediately from Lemma 17.11 Proposition 110.11 and 
Corollary O □ 

11. Directional Monotonicity II 



It is now possible to extend the directional monotonicity result of Section [S] to a 
directional monotonicity result with respect to time-space variables. 
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Proposition 11.1. Let < Amin < infQj(o) niin(l, A+, A_), h as in i6.1\} . let 

£ G (0, 1) and let f and a be the constants of Corollary \10.2[ Then each solution u 
of (GUP *^ Qi(0) such that 

J. , / i,\ ^ r Amin£^2~2 

distvvi.=(Qi(o))(w,/i) < d := -^"^ ^ 

and 

sup max(|VA+|, |9tA+|, |VA„|, \dt\A) < 6 
Qi(o) 

satisfies £~^adtu + £~^deU — |u| > in (5i/2(0) for every a G [—1, 1] and every 
e € 9i3i(0) such that ei > e; here ei denotes the first component of the vector e. 

Proof. First note that Qi H {u = 0} is by the assumptions contained in the strip 
\xi\ < of 12, implying by Corollary 110.21 and Lemma [7. II that dtu is continuous in 
Q\. We know that e~^adth + e~^deh — \h\ > 0. It follows that 

(11.1) £-^adtu + e-^deu - |u| > -SSe'^ 

provided that dist-[yi,oa(^Q^i^Qy^{u,h) < S. Suppose now towards a contradiction 
that the statement is not true. Then there exist A+,A_ G {Xmin, +oo), (t* , x*) G 
Qi/2(0)i Q!*, e*, and a solution u of (jl.l|) in Qi{0) such that dist^yi,oo(Qi(o))(u, h) < 
S, 

sup max(|VA+|, l^tA+l, |VA_|, \dt\-\) < S, 
Qi(o) 

|a*| < l,et > £ and e-^a*dtu{t* ,x*) + £~^de-u{t* ,x*) - \u{t*,x*)\ < 0. For the 
positive constant c to be defined later the functions v := £~^a*dtu + £~^d(^*u — \u\ 
and w e^^a*dtu + e~^de*u — |u| + c\x — x*]"^ — c{t — t*) satisfy then by the 
definition of 5 the following: in the set D := Qi(0) n {i; < 0} n {i < i*}, 

Au; - dtw <2nc + c- A+X{„>o} - X-X{u<o} 
+e-\X+ + X-)^^ ■ e*H"[({M = 0} n {Vm 7^ 0}) 
+e-i(A+ + X-)vta*rr [{{u = 0} n {Vu ^ 0}) 
+e"^(X{«>o}(a*9t + de')X+ - X{u<Q}{ot* dt + 5e.)A-) 
where v — (^*"^^") ^ As 

V- (a*,e*) < on {w = 0} n {u < 0} {w 0} n {e-^a*dtu + e-^ae^ii < 0} , 

we obtain by the definition of b that w is supercaloric in D provided that c has 
been chosen accordingly, say c := Ai„in/(4n). It follows that the negative infimum 
of w is attained on 

5pari? C (9pargi(0) H {i < i*}) U (Qi(0) H d{v < 0}) . 

Consequently it is attained on {t < t*} n c'paiQi(O); say at the point {t,x) G {i < 
t*} n (9par(3i(0). Since pardist((f, x), {t*,x*)) > 1/2, we obtain that 

e^^a*dtu{t,x) + e~^de'u{t,x) — \u{i,x)\ 
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= v{t, x) = W{t, x) - C\X* - iTp + c(t - t*) < -c/4 = -Aniin/(16n) . 

But this contradicts in view oi 5 ~ '^Is'"/ ^^^^- 

12. Proof of the main theorem 

The theorem is proven in several simple steps, using mainly Proposition 111.11 
and Lemma 19.11 Note that the proof can be simplified substantially in the case 
that we are dealing not with a whole class of solutions but a single solution. 
Part I: In this first part we prove uniform Lipschitz regularity and continuous 
differentiability with respect to the space variables. 

Step 1 (Directional monotonicity): Given e > 0, there are > Q and > 
(depending only on the parameters of the statement) such that 2a£~^r1dtu + 
2e~^redeU — |w| > in Qr^f-iiy) for every a e [—1,1]. The inequality holds for 
every (s,?/) e Qi/2(0) satisfying u{s,y) = 0, |VM(.s,y)| < a^re, pardist((s, y), {u > 
0}) < CTgrg and pardist((s, y), {w < 0}) < tTere, for some unit vector v^{s,y) and 
for every e G dBi satisfying e • z^e(s, ?/) > |. In particular, for e = 1, the solution u 
is by condition \1.2\ with a = (TiTi non-decreasing in Q,j/2(0) in direction (7'i,e) 
for every e e dBi{Q) such that e ■ v^{Q)> ^. 

Proof: By Lemma [9.11 there are > and > as above such that the scaled 
function Ur^{t,x) — u{s + r'^t,y + r^x)/r1 is 5 :— e^^f^fi^-close in C^((5i(0)) 
to a rotated version ft, of in Qi. Let ^^{s,!/) be the accordingly rotated version 
of the unit vector ei. Since Ur^ solves (jl.ip with respect to X+{t^ ■ +s,re ■ +y) 
and A_(r^ • +s,re ■ +y), and since max(|V(A+(r^ • +s,re • +y))\, |V(A_(7-^ • +s,re ■ 
+y))\Adt{X+{r^ ■ +s,r^ ■ +y))\,\dt{X-{rj ■ +s,r^ ■ +y))\) < Cir^, we may choose 
< 6/Ci in order to apply Proposition lll.ll to Ur^ in Qi and to conclude that 
2ae~^dtUj., + 2e~^deUr^ — > in Qi/2(0) for every a e [—1,1] and every 
e e dBi{Q) such that e • v^{s,y) > e/2. Scaling back we obtain the statement of 
Step 1. 

Step 2 (Lipschitz continuity): d{u > 0}nQri/2(0) and d{u < 0}nQri/2(0) are 
Lipschitz graphs in the direction of (0, i^ei^)) with spatial Lipschitz norms less than 

1 and temporal Lipschitz norms less than r^^ . Moreover, for each e G (0, 1) and 
{s,y) e {w = 0} n Qi/2 satisfying |Vu(s,y)| < aer., pardist((s, y), {u > 0}) < a^re 
and pardist((s, y),{u < 0}) < a^r^, the free boundaries d{u > 0} n Qr^/aC^, y) and 
d{u < 0} n Qr^/2is, y) are Lipschitz graphs (in the direction of vds, y)) with spatial 
Lipschitz norms not greater than e. 

Proof: This follows from the monotonicity obtained in Step 1. 

Step 3 (Existence of a spatial tangent plane at points (s, y) £ d{u > 

0} n d{u < 0} n (5i/2(0) satisfying |Vu(s, y)| = 0): The Lipschitz graphs of Step 

2 are both differentiable with respect to the space variables at the point (s, y), and 
the two spatial tangent planes at (s, y) coincide. 

Proof: This follows from Step 2 by letting e tend to zero. 
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Step 4 (One-phase points are regular): If {s,y) E Qri/2(0) is a free boundary 
point and the solution u is non-negative or non-positive in Qs{s,y), then the free 
boundary is the graph of a C^'"-function in Qasis, y), where ci and the C^'"-norni 
depend only on the parameters in the statement. Consequently, in Qri/2(0), there 
exist no singular one-phase free boundary points. 

Proof: By Step 2, the sets {u > 0} fl Qrj2{^) and {m < 0} n Qri/2(0) are 
sub/supergraphs of Lipschitz continuous functions. Therefore {u = 0} n Qs{s,y) 
satisfies the thickness condition required for [2l Theorem 15.1] and the statement 
follows. 

Step 5 (Existence of space normals in Qri/2{0))'- d{u > 0} n Qri/2(0) and 
d{u < 0} n (3ri/2(0) are graphs of Lipschitz continuous functions which are differ- 
entiable with respect to the space variables. 

Proof: Let {s,y) G Qri/2{0) be a free boundary point. We have to prove existence 
of a tangent plane at (s, y). 

First, if {s,y) is a one-phase point, i.e. if the solution u is non-negative or non- 
positive in Qs{s, y), then the statement holds at (s, y) by the result of Step 4. Sec- 
ond, if \Wu{s,y)\ ^ 0, the statement holds by Lemma [7.11 Last, if \Wu{s,y)\ = 
and {s,y) is the limit point of both phases {u > 0} and {u < 0}, then Step 3 
applies. 

Step 6 (Equicontinuity of the space normals): It remains to prove that the 
space normals are equicontinuous on Qri/2(0)n9{u > 0} and on Qri/2(0)n9{u < 0} 
for u in the class of solutions specified in the statement of the main theorem. 
Proof: By Step 2 we know already that the spatial Lipschitz norms of d{u > 
0} n (3,1/2(0) and d{u < 0} n Qri/2(0) are less than 1. We prove that the space 
normals are equicontinuous on (3, ^/2(0) n d{u > 0}. 

We may assume that J^(0) points in the direction of the xi-axis and that xi = 
f{t,X2,---,Xn) is the representation of d{u > 0} n Qri/2(0). Besides we have 
|V/(t,a;')| < 1 for {t,x) = {t,xi,x') e d{u > 0} n Qr,/2{0)- We claim that for 
e > there is (S^ > depending only on the parameters in the statement such that 
for any pair of free boundary points (s^, y^), (s^, y^) G d{u > 0} fl Qri/2(0), 

(12.1) pardist((si,yi),(s2,2/2)) ^ \„{s\y') - ,yi.s^y^)\ < 2e. 

In what follows let := a^r^j^ < ri/2. 

Suppose first that u is non-negative in Qp^{s^ ,y^)- Here we may as in Step 4 apply 
[2j Theorem 15.1] to the scaled function ■w{t,x) := u{s^ + p^t,y^ + p^x)/p1; since 
the C^'"-norm of the free boundary normal of w is on Qca H d{w > 0} bounded by 
a constant C3, where C2 > and C3 < -l-oo depend only on the parameters in the 
statement, we may choose 

j_ 

5e ■■= min( — ,C2)Pe 
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to obtain p27T|l . 

Next, suppose that u changes its sign at Qp, (s^,y^). If there is a point (s,y) e 
Qp^{s^ ,y^) n d{u > 0} such that \Wu{s,y)\ < then we are in the situation of 
Step 1. By Step 2 the free boundary d{u > 0}r]Qr^/2{s, y) is Lipschitz with spatial 
Lipschitz norm not greater than e. Hence ()12.1|) foUows in this case with 5^ := re/2. 
Last, if |Vu(s,y)| > pe for all points (s,y) G Qpe(s^,2/^) n 9{m > 0}, we proceed as 
follows: from the equation u{t, f{t, x'), x') = we infer that V'u + diu V'/ = on 
d{u > 0} n Qri/2(0)- Hence we obtain 



\Vf{s\{y')')^Vf{s\{y^)')\^ 



< 



(9iii(si, yi) 
\V'uis\y^)-V'u{s\y^)\ 



diu{s'^, y"^ 



Vu{s'^,y'^ 



diu{s'^,y'^) 



\diu{s\y^)\ 

\diu{s^,y'^) - 9iM(s\yi)| 



where M 



\diu{s\y^)\ 
< 4Mp-'p&idist{{s\y'),{s^,y^)) , 
I Vu|l^i/2,i(Q^^2(o)) ■ particular wc may choose 

Se 

to arrive at p2.ip . 



AM' 



Part II: We conclude the proof of the main theorem by pointing out a counter- 
example to -regularity. 

Consider the one-phase counter-example u : [— r^, r^] x [0, r] [0, +oo) from [1, p. 
376] satisfying the following: u{t,0) = for — < t < r^, and the free boundary 
touches the lateral boundary at the origin in a non-tangential way. Thus we may 
reflect u to a solution 



v(t, x) 



u{t, x), a; > 
—u{t, —x), a; < 



and obtain that is a solution of our two-phase problem in Qr for A+ = A_ = 
1. As the free boundary d{v > 0} is only Lipschitz at the origin, we conclude that 
differentiability with respect to the time variable is in general not true. □ 
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